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ABSTRACT 


In Japan, the wild house mouse Mus musculus consists of two lineages, one from Southeast Asia (Mus musculus 
castaneus, CAS) and one from northern Eurasia (Mus musculus musculus; MUS). However, the exact origins 
of the parental lineages are unclear. A recent work using mitochondrial sequences revealed that Japanese CAS 
and MUS are closely related to haplotypes from South China and the Korean Peninsula, respectively. Recent 
phylogeographic analyses using nuclear gene sequences have also confirmed a close relationship between Japan 
and Korea in the MUS component. However, the Japanese CAS components in the nuclear genome are likely 
to be unique and to differ from those of other CAS territories, including South China. Although the origins are 
still unresolved, these results allow us to conclude that two areas of the continent, South China and the Korean 
Peninsula, are the primary source areas of Japanese wild mice and suggest pre-historical introductions associated 
with certain historical agricultural developments in East Asia. 
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INTRODUCTION 


The wild house mouse, Mus musculus, occurring in the 
Japanese Islands has a short body and white belly, an outward 
appearance that differs substantially from those in the northern 
Eurasia lineage, Mus musculus musculus (MUS). Past mole- 
cular analyses have indicated that the Japanese wild mouse 
is a hybrid of the two northern and southern subspecies 
groups of MUS and M. m. castaneus (CAS) (e.g., Suzuki et 
al., 1986; Yonekawa et al., 1988, 2012). However, the exact 
geographic origins of Japanese wild mice are uncertain due 
to the veiled histories of subspecies groups in East Asia. 
From our recent phylogeographic work on mice collected 
over wide areas of the two subspecies groups (Nunome et 
al., 2010; Suzuki and Aplin, 2012; Suzuki et al., 2013), it is 
now possible to more accurately pinpoint the multiple places 
of origin of Japanese mice in the Korean Peninsula and South 
China. Here we describe recent achievements related to this 


issue and discuss the possible interpretation of the historical 
migrations of wild mice from continental areas to Japan. 


POSSIBLE SOURCE AREAS OF THE TWO 
JAPANESE LINEAGES 


Phylogeographic studies using mitochondrial cytochrome b 
(Cytb) (Fig. 1A) and control region sequences have revealed 
that M. musculus possesses five distinct subspecies lineages, 
comprising those occurring in West Europe (M. musculus 
domesticus; DOM), Nepal (un-assigned subspecies; NEP), 
and Yemen (M. musculus gentilulus), in addition to CAS 
and MUS (Bonhomme and Searle, 2012; Yonekawa et al., 
2012; Suzuki et al., 2013). CAS consists of four sublineages 
(CAS-1 to -4) that originated in the late Middle to Late Plei- 
stocene (100,000—200,000 years ago), with a trend of confined 
distribution ranges around the eastern part of the Middle East, 
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Fig. 1. Evolutionary history of mice in East Asia. (A) Cytb phylogenetic tree determined using a Bayesian MCMCMC method with the 
program BEAST v. 1.7.4. Triangles represent collapsed subclades of wild mice subspecies groups. (B) Hypothetical migration routes 
of two MUS-1 and CAS-1 sublineages. Times on the map indicate estimates for population expansions based on Cytb sequences, 
assuming a mutation rate of 8.5% per million years per lineage. Arrows on the map indicate temporal migration routes of Mus 
musculus musculus (MUS) and Mus musculus castaneus (CAS) inferred from mtDNA data. The numbers from one to six are locality 
numbers where representative haplotypes of nuclear genes were obtained as shown in Fig. 1D. (C) A schematic illustration for intact 
haplotype degradation of MUS (blue bar) and CAS (yellow bar) through generations by chromosomal recombination, assumed to 
occur at a rate of 0.52 cM/Mb (Jensen-Seaman et al., 2004). This rate means one recombination event in a 200 kb genomic region 
in 1,040 generations. Thus, it can be assumed that 5,200 generations will produce five recombinations in a 200 kb genomic region. 
(D) Representative haplotypes from China, Korea, and Japan discovered in Nunome et al. (2010). The letters “m”, “c”, and “d” in 
boxes represent MUS-, CAS-, and Mus musculus domesticus (DOM)-type alleles of eight linked genes within a 200 kb region on 
mouse chromosome 8. The letter “u” refers to the alleles that could not be assigned to a particular type of subspecies. Haplotypes 


in Japan contain genomic fragments from three subspecies that suggest ancient contacts of the subspecies there. 


except for CAS-1, which appears to have spread rapidly 
during prehistoric times from a source region somewhere in 
India to the far eastern periphery of the CAS territory, includ- 
ing South China and Japan (Suzuki et al., 2013). Notably, 
the CAS-1 sublineage has low nucleotide diversity, which is 
consistent with the recent phylogeographic work on the 
castaneus subspecies group by Rajabi-Maham et al. (2012); 
see also Bonhomme and Searle (2012). Notably, network 
analyses depict a sub-cluster (“CAS-1a”) that harbors haplo- 
types from northern Japan (Hokkaido and northern Honshu), 
South China (Kunming and Guilin), and South Sakhalin. The 
presence of the subgroups is suggestive of a stepwise histori- 
cal range expansion of CAS-1. The expansion time of the 
Cytb haplotypes, at a relatively higher mutation rate (e.g., 
8.5% per million years per lineage), was roughly estimated 
to be 9,000 and 4,500 years ago for CAS-1 and CAS-la, 
respectively (Suzuki et al., 2013; Kuwayama et al., unpublish- 
ed) (Fig. 1B). 
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We suspect that the expansion of CAS-1 mice reflected 
ecological transformation of the landscape by early agricultur- 
alists and the emergence of urban centers and trade networks, 
as in the case of the Middle East populations of M. musculus 
(Auffray et al., 1990; Cucchi and Vigne, 2006). There is 
good evidence for domestication of cereal crops, including 
rice and millet, by about 9,000 years ago in South and East 
Asia (Khush, 1997; Londo et al., 2006; Liu et al., 2007; 
Molina et al., 2011). Thus, the CAS-1 range expansion may 
be associated with the spread of cereal cultivation in these 
regions (Suzuki et al., 2013). With regard to the CAS-la 
group expansion, it is of interest to note recent molecular 
and archaeological evidence for the origin of cultivated rice 
(Fuller et al., 2010, 2011; Huang et al., 2012). Archaeological 
studies have found that 4,000 years ago rice cultivation had 
spread to the nearby Yunnan province where Kunming is 
located (Fuller et al., 2010, 2011). An intensive genome sur- 
vey on wild and cultivated rice suggested that southern China 
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is the place of origin of cultivated rice in Asia (Huang et al., 
2012). Thus, it is reasonable to believe that the CAS-la group 
expanded its range in association with the spread of rice cul- 
ture from South China to a wide area of East Asia including 
the Japanese Islands. 

The MUS lineage consists of two sublineages, major (MUS- 
1) and minor (MUS-2) groups, which cover a broad range of 
North China and Russia, with an estimated divergence time 
of 150,000 years ago (Suzuki et al., 2013). Cytb sequence 
analyses have provided evidence for rapid population expan- 
sion of MUS-1, with an estimated time of 26,000 years ago 
at a mutation rate of 8.5% per million years per lineage 
(Suzuki et al., 2013) (Fig. 1B). This suggests an ancient resi- 
dence of the MUS mice in north China and the Korean Penin- 
sula, before the predicted ancient expansion of CAS-1 in the 
northern part of its range. In addition, the MUS-1 sublineage 
can be further subdivided into three subgroups. One subgroup, 
“MUS-Ic,” exclusively integrated haplotypes from Japan, 
Korea and its nearby continental areas, northeast China, and 
southern Primorye, Russia. It is interesting to note that the 
estimated expansion time of MUS-Ic, about 7,700 years 
ago, is somewhat comparable to that of CAS-1 (Suzuki et 
al., 2013). 

Our mitochondrial DNA analyses clearly indicated that 
the mitochondrial DNA of Japanese wild mice originated 
from southern China and the Korean peninsula, for CAS and 
MUS, respectively. However, we must carefully consider the 
timing of the introduction of each of these two components 
into Japan. To investigate the order of the migration of Japan- 
ese MUS and CAS lineages, we calculated the mean genetic 
distances between haplotypes from the predicted source areas 
of the continent and Japan using the Cytb sequences with 
Kimura-2-parameter (K2P) genetic distance. The mean gene- 
tic distances between the continental and Japanese groups 
were 0.02 and 0.001 for the CAS and MUS subgroups, res- 
pectively (Kuwayama et al., unpublished). Thus, we suggest 
that CAS arrived in Japan prior to MUS. The recent MUS 
migration from the Korean peninsula to Japan is likely to be 
associated with the peopling of Japan a few thousand years 
ago. Phylogenetic analyses using Y chromosomal and mito- 
chondrial markers of Japanese people have consistently indi- 
cated an immigration of Japanese people from Korea during 
the Neolithic period (Tanaka et al., 2004; Hammer et al., 
2006). 


NUCLEAR GENE COMPONENT SPECIFICITY 
IN JAPANESE WILD MICE 


Our previous phylogeographic study focusing on haplotypes 
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of linked autosomal genes highlighted the concealed evolu- 
tionary history of Japanese wild mice (Nunome et al., 2010). 
In that study, we examined eight linked genes in a 200 kb 
region of chromosome 8 and assigned the alleles into three 
subspecies groups (DOM, CAS, and MUS) and disclosed 
haplotype structures defined by the combination of the eight 
genes and their geographic distributions. The majority of 
mice from Japan have a MUS-specific haplotype. Interesting- 
ly, a minor proportion of mice, particularly those from north- 
ern Japan, have “recombinant” haplotypes, in which CAS 
alleles or CAS+DOM alleles are embedded in the 200-kb 
MUS haplotypes. 

The existence of small CAS segments within the MUS- 
like haplotypes implies that a number of generations existed 
after the hybridization events between the CAS and MUS 
groups in Japan (Fig. 1C). Hybridization time has been esti- 
mated to be a thousand generations (Nunome et al., 2010), 
providing concrete evidence that the “recombinant” haplo- 
types resulted from a historical event in ancient times and 
were not due to stowaway introduction during the modern 
age. The results of mitochondrial and nuclear gene analyses 
allowed us to predict that two areas of the continent, South 
China and the Korean peninsula, are the primary source areas 
of Japanese wild mice. A CAS subgroup would have arrived 
in Japan from southern China as a precursor to mice in Japan. 
Then, a MUS subgroup subsequently entered Japan from the 
Korean peninsula, perhaps a few thousand years ago, and 
mixed with the existing Japanese CAS populations. 

Notably, the Japanese CAS components in the nuclear 
genome are likely to be unique and to differ from those of 
other CAS territories, including South China. Allele types of 
Dbndd1, for example, differ between northern Japan (allele 
“c2”) and other areas (allele “c1”) (Fig. 1D). Phylogeographic 
analyses using nuclear genes have revealed the presence of 
the East European DOM within a CAS segment of the haplo- 
types from northern Japan, implying the presence of hidden 
evolutionary episodes of the CAS migration from the conti- 
nent to Japan (Nunome et al., 2010) (Fig. 1D). Thus, to fur- 
ther understand the histories of mice in eastern Eurasia, an 
intensive geographic survey to identify the unassigned CAS 
and DOM fragments incorporated into the haplotypes from 
northern Japan is desirable. 

The ranges of the three M. musculus subspecies, particu- 
larly in the eastern part of Eurasia, have not been finely deter- 
mined and their evolutionary history remains unclear. Phy- 
logeographic analyses of wild mice from the Japanese Islands, 
the terminal regions of several different migrations from 
west to east, and those from the two source areas of South 
China and the Korean Peninsula, provide a valuable founda- 
tion for building a comprehensive view on the evolutionary 
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history of the house mouse, Mus musculus, and highlight 
hidden episodes of ancient peopling and agricultural develop- 
ment in East Asia as well. 


ACKNOWLEDGMENTS 


We thank Kuniya Abe, Ken Aplin, Sang-Hoon Han, Naoto 
Hanzawa, Hidetoshi Ikeda, Mei-Lei Jin, Gohta Kinoshita, 
Alexey P. Kryukov, Yoichi Matsuda, Yoshifumi Matsushima, 
Ibnu Maryanto, Nobumoto Miyashita, Kimiyuki Tsuchiya, 
Peter Vogel, Toshihiko Shiroishi, Toyoyuki Takada, and 
Hiromichi Yonekawa for their valuable comments, and 
numerous other mouse collectors for kind help in supplying 
the valuable samples used in this study. This study was sup- 
ported in part by a Grant-in-Aid for Scientific Research (B) 
from Japan Society for the Promotion of Science (JSPS, 
2440513). 


REFERENCES 


Auffray JC, Vanlerberghe F, Britton-Davidian J, 1990. The house 
mouse progression in Eurasia: a palaeontological and arch- 
aeozoological approach. Biological Journal of the Linnean 
Society, 41:13-25. 

Bonhomme F, Searle JB, 2012. House mouse phylogeography. 
In: Evolution of the house mouse, Cambridge series in 
morphology and molecules (Eds., Macholan M, Baird SJE, 
Munclinger P, Pidlek J). Cambridge University Press, Cam- 
bridge, pp. 278-296. 

Cucchi T, Vigne JD, 2006. Origin and diffusion of the house 
mouse in the Mediterranean. Human Evolution, 21:95-106. 

Fuller DQ, Sato YI, Castillo C, Qin L, Weisskopf AR, Kingwell- 
Banham EJ, Song J, Ahn SM, van Etten J, 2010. Consilience 
of genetics and archaeobotany in the entangled history of 
rice. Archaeological and Anthropological Sciences, 2:115- 
131. 

Fuller DQ, van Etten J, Manning K, Castillo C, Kingwell-Banham 
E, Weisskopf A, Qin L, Sato YI, Hijmans RJ, 2011. The 
contribution of rice agriculture and livestock pastoralism to 
prehistoric methane levels: an archaeological assessment. 
The Holocene, 21:743-759. 

Hammer MF, Karafet TM, Park H, Omoto K, Harihara S, Stonek- 
ing M, Horai S, 2006. Dual origins of the Japanese: common 
ground for hunter-gatherer and farmer Y chromosomes. 
Journal of Human Genetics, 51:47-58. 

Huang X, Kurata N, Wei X, Wang ZX, Wang A, Zhao Q, Zhao 
Y, Liu K, Lu H, Li W, Guo Y, Lu Y, Zhou C, Fan D, Weng 
Q, Zhu C, Huang T, Zhang L, Wang Y, Feng L, Furuumi 
H, Kubo T, Miyabayashi T, Yuan X, Xu Q, Dong G, Zhan 
Q, Li C, Fujiyama A, Toyoda A, Lu T, Feng Q, Qian Q, Li 
J, Han B, 2012. A map of rice genome variation reveals the 


270 


origin of cultivated rice. Nature, 490:497-501. 

Jensen-Seaman MI, Furey TS, Payseur BA, Lu Y, Roskin KM, 
Chen CF, Thomas MA, Haussler D, Jacob HJ, 2004. Com- 
parative recombination rates in the rat, mouse, and human 
genomes. Genome Research, 14:528-538. 

Khush GS, 1997. Origin, dispersal, cultivation and variation of 
rice. Plant Molecular Biology, 35:25-34. 

Liu L, Lee GA, Jiang L, Zhang J, 2007. Evidence for the early 
beginning (c. 9000 cal. BP) of rice domestication in China: 
a response. The Holocene, 17:1059-1068. 

Londo JP, Chiang YC, Hung KH, Chiang TY, Schaal BA, 2006. 
Phylogeography of Asian wild rice, Oryza rufipogon, reveals 
multiple independent domestications of cultivated rice, Oryza 
sativa. Proceedings of the National Academy of Sciences 
of the United States of America, 103:9578-9583. 

Molina J, Sikora M, Garud N, Flowers JM, Rubinsteina S, 
Reynolds A, Huang P, Jackson S, Schaal BA, Bustamante 
CD, Boyko AR, Purugganan MD, 2011. Molecular evidence 
for a single evolutionary origin of domesticated rice. Proceed- 
ings of the National Academy of Sciences of the United 
States of America, 108:8351-8356. 

Nunome M, Ishimori C, Aplin KP, Tsuchiya K, Yonekawa H, 
Moriwaki K, Suzuki H, 2010. Detection of recombinant hap- 
lotypes in wild mice (Mus musculus) provides new insights 
into the origin of Japanese mice. Molecular Ecology, 19: 
2474-2489. 

Rajabi-Maham H, Orth A, Siahsarvie R, Boursot P, Darvish J, 
Bonhomme F, 2012. The south-eastern house mouse Mus 
musculus castaneus (Rodentia: Muridae) is a polytypic sub- 
species. Biological Journal of the Linnean Society, 107:295- 
306. 

Suzuki H, Aplin KP, 2012. Phylogeny and biogeography of the 
genus Mus in Eurasia. In: Evolution of the house mouse 
(Eds., Macholan M, Baird SJE, Munclinger P, Pidlek J). 
Cambridge University Press, Cambridge, pp. 35-64. 

Suzuki H, Miyashita N, Moriwaki K, Kominami R, Muramatsu 
M, Kanehisa T, Bonhomme F, Petras ML, Yu ZC, Lu DY, 
1986. Evolutionary implication of heterogeneity of the non- 
transcribed spacer region of ribosomal DNA repeating units 
in various subspecies of Mus musculus. Molecular Biology 
and Evolution, 3:126-137. 

Suzuki H, Nunome M, Kinoshita G, Aplin KP, Vogel P, Kryukov 
AP, Jin ML, Han SH, Maryanto I, Tsuchiya K, Ikeda H, 
Shiroishi T, Yonekawa H, Moriwaki K, 2013. Evolutionary 
and dispersal history of Eurasian house mice Mus musculus 
clarified by more extensive geographic sampling of mito- 
chondrial DNA. Heredity, advance online publication. http:// 
dx.doi.org/10.1038/hdy.2013.60. 

Tanaka M, Cabrera VM, González AM, Larruga JM, Takeyasu 
T, Fuku N, Guo LJ, Hirose R, Fujita Y, Kurata M, Shinoda 
K, Umetsu K, Yamada Y, Oshida Y, Sato Y, Hattori N, 
Mizuno Y, Arai Y, Hirose N, Ohta S, Ogawa O, Tanaka Y, 
Kawamori R, Shamoto-Nagai M, Maruyama W, Shimokata 
H, Suzuki R, Shimodaira H, 2004. Mitochondrial genome 
variation in Eastern Asia and the peopling of Japan. Genome 


Anim. Syst. Evol. Divers. 29(4), 267-271 


Research, 14:1832-1850. 

Yonekawa H, Moriwaki K, Gotoh O, Miyashita N, Matsushima 
N, Shi LM, Cho WS, Zhen XL, Tagashira Y, 1988. Hybrid 
origin of Japanese mice “Mus musculus molossinus”: evid- 
ence from restriction analysis of mitochondrial DNA. Mole- 
cular Biology and Evolution, 5:63-78. 

Yonekawa H, Sato JJ, Suzuki H, Moriwaki K, 2012. Origin 
and genetic status of Mus musculus molossinus: a typical 


Anim. Syst. Evol. Divers. 29(4), 267-271 


Evolutionary History of the Japanese House Mouse 


example for reticulate evolution in the genus Mus. In: Evolu- 
tion of the house mouse (Eds., Macholan M, Baird SJE, 
Munclinger P, Pidlek J). Cambridge University Press, Cam- 
bridge, pp. 94-113. 


Received April 9, 2013 
Revised August 26, 2013 
Accepted August 31, 2013 


271 


